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Abstract. The 6.4 s X-ray pulsar IE 1048.1-5937 was ob- 
served by BeppoSAX in 1997 May. This source belongs to 
the class of "anomalous" pulsars which have pulse periods 
in range 5-11 s, show no evidence of optical or radio coun- 
terparts, and exhibit long-term increases in pulse period. 
The phase-averaged 0.5-10 keV spectrum can be described 
by an absorbed power-law and blackbody model. The 
best-fit photon index is 2.5±0.2 and the blackbody tem- 
perature and radius are 0.64±0.01 keV and 0.59 ±0.02 km 
(for a distance of 3 kpc), respectively. The detection of 
blackbody emission from this source strengthens the simi- 
larity with two of the more well studied "anomalous" pul- 
sars, IE 2259+586 and 4U 0142+614. There is no evidence 
for any phase dependent spectral changes. 

The pulse period of 6.45026 ± 0.00001 s implies that 
IE 1048.1-5937 continues to spin-down, but at a slower 
rate than obtained from the previous measurements in 
1994 and 1996. 



pulsars: individual 
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1. Introduction 



IE 1048.1-5937 is an unusual pulsar, with a period of 6.4 s 
and a soft spectrum, discovered during Einstein obser- 
vations of the Carina nebula (Seward et al. 1983). The 
source has been spinning down for at least the last 17 
years (Mereghetti 1995; Corbet & Mihara 1997). The spec- 
trum has been modeled by an absorbed power-law with a 
photon index, a, of ~2-3 (Seward et al. 1983; Corbet & 
Mihara 1997). This spectral shape is softer than those of 
typical high-luminosity X-ray pulsars. Despite a small er- 
ror box, no optical counterpart has been identified, with a 
limiting magnitude of my ~ 20 which excludes the pres- 
ence of a massive companion (Mereghetti et al. 1992). A 
recent observation with the RossiXTE satellite provides a 
strong upper limit to the projected X-ray semi-major axis 



of 0.06 lt-s for orbital periods between 200 s and ~1 day 
(Mereghetti et al. 1997). The lack of an optical counterpart 
and orbital Doppler shifts argue against a binary model 
for IE 1048.1-5937, unless the companion has a low mass. 
Mereghetti et al. (1997) find that a probable upper limit 
to the mass of a Roche lobe filling main sequence com- 
panion is ~0.3 M . Masses up to ~0.8 M Q are allowed in 
the case of a helium-burning companion filling its Roche 
lobe. 

The properties described above are similar to those of 
a small number of other X-ray pulsars with spin periods 
in the 5-11 s range (Mereghetti & Stella 1995), such as 
IE 2259+586 and 4U 0142+614. These form a class of so- 
called "anomalous" pulsars, with clearly different proper- 
ties from the majority of systems. Although accretion from 
a very low mass companion cannot be excluded, the lack 
of evidence for a binary nature from any of these systems 
has stimulated models where the X-ray emission originates 
from a compact object that is not in an interacting binary 
system. While an isolated, massive, white dwarf powered 
by the loss of rotational energy, as originally proposed for 
IE 2259+586 (Paczynski 1990; Usov 1994), has been ruled 
out by the detection of a large increase in the spin-down 
rate of IE 1048.1-5937 (Mereghetti 1995), other single ob- 
ject models, such as loss of magnetic energy of a strongly 
magnetized neutron star (Thompson & Duncan 1993), or 
an isolated neutron star accreting from a circumstcllar 
disk (Corbet et al. 1995; van Paradijs et al. 1995), may be 
applicable. 

We present a study of IE 1048.1-5937 based on data 
obtained with the BeppoSAX satellite. We focus on the 
X-ray spectrum at energies <10 keV and on the pulse 
period history. As with some of the other "anomalous" 
pulsars, we find evidence for the presence of a blackbody 
spectral component. Since this component is not observed 
from the majority of other accreting X-ray pulsars, this 
strengthens the similarity between IE 1048.1-5937 and the 
other better studied "anomalous" pulsars. 
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2. Observations 

Results were obtained with the Low-Energy Concentrator 
Spectrometer (LECS; 0.1-10 keV; Parmar et al. 1997a) 
and Medium-Energy Concentrator Spectrometer (MECS; 
1.3-10 keV; Boella et al. 1997) instruments. The MECS 
consists of three identical grazing incidence telescopes 
with imaging gas scintillation proportional counters in 
their focal planes. The LECS uses an identical telescope 
as the MECS, but utilizes an ultra-thin (1.25 /xm) en- 
trance window and a driftless configuration to extend the 
low-energy response to 0.1 keV. The fields of view of the 
LECS and MECS are circular with diameters of 37' and 
56', respectively. In the overlapping energy range, the po- 
sition resolution of both instruments is similar and corre- 
sponds to 90% encircled energy within a radius of 2'5 at 
1.5 keV. At lower energies, the encircled energy is propor- 
tional to E" - 5 . In the central 10' the LECS 0.1-10 keV 
and the MECS 1-11 keV background counting rates are 
9.7 x 10~ 5 arcmin" 2 s" 1 and 10.9 x 10" 5 arcmin" 2 s" 1 
(for 2 MECS units), respectively. 

IE 1048.1-5937 was observed by BeppoSAX on 1997 
May 10 01 h 24 m to May 11 15 h 24 m UTC. Good data were 
selected from time intervals when the minimum elevation 
angle above the Earth's limb was >4° and when the in- 
strument configurations were nominal using the SAXDAS 
1.2.0 data analysis package. This gives exposures of 80 ks 
for the MECS and 32 ks for the LECS, which was only op- 
erated during satellite night-time. One of the three MECS 
units failed on 1997 May 9 (one day before our observa- 
tion), and data were only obtained from the remaining 
two MECS units. 

3. Results 

3.1. Spectral analysis 

Spectra were obtained centered on the pulsar position us- 
ing extraction radii of 8' and 4' for the LECS and MECS, 
respectively. The appropriate LECS response matrix was 
generated using LEMAT 3.4.0. Background subtraction 
was performed using standard blank field exposures. The 
background subtracted IE 1048.1-5937 count rates are 
0.08 and 0.12 s" 1 in the LECS and MECS, respectively. 
Examination of the LECS spectrum reveals that the pul- 
sar is only detected above 0.5 keV, and data below this en- 
ergy are excluded. Similarly, the MECS fits are restricted 
to the energy range 1.65-10 keV. All the spectra were re- 
binned to have >20 counts in each bin to allow the use of 
X 2 statistics. 

In order to compare our results with previous obser- 
vations, the spectra were first fit with two models: (i) an 
absorbed power-law, and (ii) an absorbed power-law plus 
blackbody. Although a single power-law with a=4.04 sat- 
isfactorily describes the MECS spectrum, it is clearly in- 
adequate when the LECS spectrum is included, giving a 
xl of 1.520 for 304 degrees of freedom (dof). Figure |] 



Table 1. Spectral fit results for IE 1048.1-5937. Uncer- 
tainties are given at 68% confidence for one parameter of 
interest. The normalization is at 1 keV 



Parameter 


Value 


power-law model 




a 


3.36 ±0.05 


N H (10 22 atoms cm" 2 ) 


1.54 ±0.10 


xS(dof) 


1.520 (304) 


power-law & blackbody model 




Oi 


o co _i_ n on 


k i bb V ) 


n ra -t- n m 


±tbb l Km J 


n A- n no 


1\T TT /in^2 r,f A mc nm~ ^\ 

L > H iru atoms cm i 




Xv( aot ) 


f) QQ7 ( QflO^ 


Broken power-law model 




Ctl 


1.71 ± 0.2 


O 2 


o.uo in u.uo 


Break Energy (keV) 


1 71 _l_ A O 

1. 1 1 ± 0.2 


N H (10 22 atoms cm" 2 ) 


0.47 ± 0.07 




1.015 (302) 


Two blackbody model 




kTbbl (keV) 


0.58 ±0.02 


Rbbl a (km) 


0.77 ±0.08 


kT bb 2 (keV) 


1.25 ±0.12 


Rbb2 a (km) 


0.10 ±0.03 


N H (10 22 atoms cm" 2 ) 


0.14 ±0.04 


xS(dof) 


0.964 (302) 


Cut-off power-law model 




Q 


0.42 ±0.07 



Cut-off energy (keV) 1.43 ± 0.13 

N H (10 22 atoms cm" 2 ) 0.42 ± 0.07 
X 2 (dof) 1.041 (303) 

"For a distance of 3 kpc 



shows the fit residuals before and after the inclusion of 
the blackbody component. The power-law plus blackbody 
model gives a xl OI 0.997 for 302 dof for a blackbody 
temperature, KTbb, of 0.64 ± 0.01 keV and radius, Rbb, 
of 0.59 ± 0.02 km for a distance of 3 kpc. The best-fit 
value of a is 2.52 ± 0.20 and low-energy absorption, Nh, 
equivalent to (4.5 ± 1.0) x 10 21 atoms cm" 2 is required. 
An F-test shows that the inclusion of the blackbody is 
highly significant (the reduction in x 2 has a formal prob- 
ability of ^10" 28 , corresponding to ^llcr). The best fit 
two-component spectrum is plotted in Fig. [Hand the cor- 
responding parameters are given in Table TO. The unab- 
sorbed 2-10 keV source flux is 7.0xl0~ 12 ergs cm" 2 s" 1 , 
with the blackbody component accounting for 55% of the 
flux in this energy range. For comparison, the results of 
fitting the same spectral model to some other "anoma- 
lous" pulsars are given in Table ^. In order to investigate 
whether other simple spectral models can also give accept- 
able representations of the IE 1048.1-5937 spectrum, fits 
using three additional models were also performed (see 
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Table 2. The pulsed fraction as a function of energy as 
determined from the MECS data 



Energy Range 


Pulsed Fraction 


(keV) 


(background corrected) 


1.0-1.5 


0.76±0.08 


1.5-2.0 


0.74±0.05 


2.0-2.5 


0.77±0.04 


2.5-4.0 


0.84±0.03 


4.0-6.0 


0.79±0.05 


6.0-10.0 


0.54±0.12 



Table |lj). All these models gave fits of comparable quality 
as the power-law and blackbody combination. This means 
that we are unable to meaningfully distinguish between 
the different spectral models listed in Table [|, with the 
exception of the power-law. However, in order to compare 
our results with previous ones, we will use power-law and 
blackbody fits in the subsequent analysis. 




Phase 



Fig. 2. The 0.5-10 keV pulse profile of IE 1048.1-5937. 
For clarity, the pulse profile is repeated 



3.2. Pulse timing and phase resolved spectroscopy 

The MECS counts were used to determine the IE 1048.1- 
5937 pulse period, after correction of their arrival times to 
the solar system barycenter. The data were divided into 
10 time intervals (each with ^950 counts) and for each in- 
terval the relative phase of the pulsations determined. The 
phases of the 10 time intervals were then fit with a linear 
function giving a best-fit period of 6.45026 ± 0.00001 s. 
This value is plotted together with the previous measure- 
ments in Fig. |. The 0.5-10 keV folded light curve of 
IE 1048.1-5937 shows an approximately sinusoidal vari- 
ation with an amplitude of ^70% (Fig. ||). The pulsed 
fraction as a function of energy was determined in the fol- 
lowing way: The data for each energy range was folded at 
the best-fit period and fit with a model consisting of a con- 
stant, C, and a sine term (f(<j>) = C + A sm((j> + (j>o)) . The 
only free parameters in this fit are C and the amplitude 
of the sine, A, since 4>o was determined from the total en- 
ergy range data and was fixed during the fits (there is no 
evidence for a change of tfio with energy) . We then deter- 
mined the pulsed fraction denned as the amplitude of the 
sine term (A) divided by C (a correction for background 
was applied). The results are given in Table ||. The cor- 
rection for background is only important in the highest 
energy range. From this table we conclude that the pulsed 
fraction is approximately constant over the whole energy 
range. We caution that the systematic errors are largest in 
the highest energy bin (with the background contributing 
33% of the mean count rate). 

A set of 4 phase-resolved spectra of the pulsar were 
accumulated. These were then fit with the power-law plus 
blackbody model used in Sect. 3.1, with Nh fixed at the 



phase-averaged best-fit value. There are insufficient counts 
to simultaneously constrain both the power-law and black- 
body components. Initially, the blackbody spectral param- 
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Fig. 3. The pulse period of IE 1048.1-5937 as a function 
of time. Data have been obtained from Mereghetti (1995), 
Corbet & Mihara (1997) and Mereghetti et al. (1997) 



eters were fixed at their phase-averaged best-fit values and 
only the power-law parameters were allowed to vary. The 
resulting fits are unacceptable, primarily because the spec- 
tra obtained at the valleys are not well fit (x^=2.2 for 79 
dof) due to the fixed blackbody component contributing 
too much flux. The fits were then repeated with the black- 
body normalization allowed to vary, giving acceptable val- 
ues of xl (second panel of Fig. |^). 

The best-fit values of a obtained with a variable black- 
body normalization are shown in Fig. ^ (bottom panel) 
and reveal a marginal phase dependence (% 2 of 11.1 for 3 
dof, corresponding to the 1% level). A phase dependence of 
the blackbody radius is evident. The variation in a may 
be caused by inadequacies in the model, or correlations 
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power— law power— law + blackbody 




Fig. 1. The spectrum of IE 1048.1-5937 obtained with the LECS and two MECS units. The left panels show the 
spectrum and the absorbed power-law model fit. Note the large residuals below 1.5 keV. The right panels show the 
spectrum with the absorbed power-law plus blackbody model fit. The units of the residuals are counts s _1 keV -1 

Table 3. The period history for IE 1048.1-5937 



Time 


Period 


pa 


Pulsed Fraction 


Observatory 


Reference 


(JD-2440000) 


(s) 


10" 4 s yr" 1 








4068.7 


6.4377±0.001 




0.68±0.07 


Einstein 


Seward et al. (1986) 


6236.6 


6.4407±0.0009 


5±2 


0.68 


EXOSAT 


Seward et al. (1986) 


7267.0 


6.44153±0.0002 






Ginga 


Corbet & Day (1990) 


7277.0 


6.44185±0. 00001 


4±3 


b 


Ginga 




7433.0 


6.4422±0.0004 






Ginga 




8787.24 


6.444868±0.000007 


7.29±0.03 


~0.70 


ROSAT 


Mereghetti (1995) 


8844.48 


6.44499±0. 000034 






ROSAT 


!) 


8972.74 


6.44532±0. 000072 






ROSAT 


If 


8994.80 


6.445391±0.000013 






ROSAT 


If 


9417.0 


6.446645±0. 000001 


10.30±0.04 


0.70±0.05 c 


ASCA 


Corbet & Mihara (1997) 


10294.5 


6.449769±0.000004 


12.99±0.02 


d 


RossiXTE 


Mereghetti et al. (1997) 


10579.31 


6.45026±0.000013 


6.3±0.2 


0.76±0.02 


BeppoSAX 


This paper (see also Table 0) 



a The values of P are calculated by comparing the best pulse period determined by each satellite with respect 
to that of the previous one, i.e. a long term average. 

No value for the pulsed fraction is reported for the Ginga observation due to source confusion. 
c Obtained from Fig. 2 in Corbet & Mihara (1997) 
d Not reported. 



between the derived fit-parameters. The phase-resolved 
spectra were also fit with both a and kTbb fixed and both 
normalizations free. This also gives acceptable values of 
xl which are almost identical to those obtained from the 
fit described above. The derived blackbody radii are not 
significantly different from those obtained with the vari- 



able power-law index fits (see Fig. We conclude that 
there is no evidence for any phase dependent changes in 
either a or kTbb- 
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P 


kT bb 


Rbb 


Lbb 


Q 


ra 

pulse 


References 


Source 


00 


(keV) 


(km) 


Ltot 








4U 0142+61 


8.69 


0.39 


2.4di kpc 


0.4 


3.67+0.09 


-10% 


White et al. (1996) 


RX J0720.4-3125 6 


8.38 


0.079±0.004 c 










Haberl et al. (1996) 


IE 1048-5937 


6.45 


0.64+0.01 


0.59+0.02 d 3kpc 


0.55 


2.5+0.2 


-70% 




1RXS J170849.0-400910 


11.00 


0.41+0.03 


4+0.6 diokpc 


0.17 


2.9+0.3 


-30% 


Sugizaki et al. (1997) 


IE 1841-045 


11.77 








3.0+0.2 


-35% 


Vasisht & Gotthelf (1997) 


IE 2259+586 


6.98 


0.44+0.01 


3.3+0.3 d 4 k P c 


0.4 


3.93+0.09 


-30% 


Corbet et al. (1995); 
Parmar et al. (1997b) 



a Amplitude of pulsed emission 

h Membership of "anomalous" pulsar class uncertain 
c Obtained from one component blackbody fit 
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Fig. 4. Spectral fit parameters as a function of pulse 
phase. Errors are at 68% confidence for 1 parameter of 
interest. The blackbody radius is for an assumed distance 
of 3 kpc, the flux in the top-panel is in units of 10~ 12 ergs 
s -1 cm~ 2 . Pulse phase is the same as in Fig. @ 



4. Discussion 

The 0.5-10 keV spectrum of IE 1048.1-5937 can be de- 
scribed by the sum of an absorbed power-law and black- 
body models. The existence of a blackbody component in 
IE 1048.1-5937 was first suggested by the ASCA results 
of Corbet & Mihara (1997). However, these authors are 
unable to clearly discriminate between this two compo- 
nent model, which gives a \l °f 0-94 for 332 dof, and an 
absorbed power-law which also provides an acceptable de- 
scription of the ASCA data with a xl of 1.02 for 337 dof. 
The BeppoSAX results reported here clearly require that 
the IE 1048.1-5937 spectrum differs significantly from an 
absorbed power-law. This deviation in consistent with a 
blackbody, but we cannot exclude the possibility that it 
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Fig. 5. The 2-10 keV flux of IE 1048.1-5937 as a function 
of time. Data are from Seward et al. (1986), Mereghetti 
(1995) and Corbet & Mihara (1997). The ASCA and Bep- 
poSAX values are for a power-law fit (to be consistent 
with the other measurements) ; the values obtained with 2- 
component fits are —15% lower. The first point (open tri- 
angle) is the upper limit reported in Seward et al. (1986), 
where a short discussion about the uncertainties in this 
value can be found. ROSAT and Einstein fluxes in the 2- 
10 keV band are taken from Mereghetti (1995) and Seward 
et al. (1986) 

has another form. A discussion about the physical inter- 
pretation of the alternate models can be found in White 
et al. (1996). This discussion is applicable to IE 1048.1- 
5937, since its spectral shape is roughly similar to that 
of 4U 0142+614. The reason for the significant BeppoSAX 
detection of spectral complexity is probably related to the 
combination of good energy resolution and extended low 
energy coverage of the LECS, together with the longer 
BeppoSAX exposure. 
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IE 2259+586 (Corbet et al. 1995) and 4U 0142+614 
(White et al. 1996) have also been successfully fit with 
absorbed power-law and blackbody spectral models. Since 
the majority of X-ray pulsar spectra are fit by absorbed 
power-law models in the 0.5-10 keV energy range (e.g., 
White et al. 1983), our results strengthen the similarity 
between IE 1048.1-5937 and the other "anomalous" pul- 
sars. Table [| is a compilation of the spectral properties of 
these "anomalous" pulsars. The blackbody component in 
these systems has been interpreted as evidence for quasi- 
spherical accretion onto an isolated neutron star formed 
after common envelope evolution and spiral-in of a mas- 
sive X-ray binary (White et al. 1996). In this case, the 
accretion flow results from the remaining part of the mas- 
sive star's envelope and may consist of two components 
(Ghosh et al. 1997). A low-angular momentum component 
gives rise to the blackbody emission from a considerable 
fraction of the neutron star surface, while a high-angular 
momentum component forms an accretion disk and is re- 
sponsible for the power-law emission and the long term 
spin-down evolution. 

The area for the blackbody emitting surface obtained 
for IE 1048.1-5937 (-0.59 d 3kpc km) is smaller than those 
of the other "anomalous" systems. However, the distance 
to IE 1048.1-5937 is poorly constrained and the assumed 
distance of 3 kpc could be considered as a lower limit 
since the measured Nh implies that it lies behind the Ca- 
rina Nebula at 2.8 kpc (Seward et al. (1986). White et al. 
(1996) propose that the low pulsed fraction observed from 
4U 0142+614 results from the large polar cap area in this 
system. This is consistent with the small polar cap area 
and high pulsed fraction reported here for IE 1048.1-5937 
but not with IE 2259+586 and 1RXS J170849.0-400910, 
which both have large radii and a moderate pulsed frac- 
tion (see Table f|). Interestingly, the best-fit blackbody 
temperature of 0.64 keV is somewhat higher than for 
IE 2259+586 and 4U 0142+614 which may be related to 
the smaller area of the polar caps. Summarizing, we find 
that for IE 1048.1-5937 (i) the blackbody temperature is 
higher, (ii) the blackbody radius is smaller, (iii) the power- 
law index is smaller (i.e. the spectrum is harder), and (iv) 
the pulsed fraction is higher than for the other "anoma- 
lous" pulsars (with the exception of RX J0720. 4-3125 
where the blackbody parameters are obtained from a one 
component fit). 

The phase-resolved spectra are not well fit with a con- 
stant contribution from the blackbody component. This 
is unsurprising given the large pulse amplitude (—70%), 
together with the large blackbody contribution (55%; 2- 
10 keV) to the phase averaged flux. The probable con- 
stancy of the pulsed fraction over the 1-10 keV energy 
range and the lack of any spectral dependence on pulse 
phase (see Sect. 3.2) imply that the whole spectrum is 
varying in a similar manner; i.e. the pulsed component 
cannot be attributed solely to either the blackbody or the 
power-law components. Furthermore, the long term vari- 



ations in source flux (Fig. ||) and the approximate con- 
stancy of the pulsed fraction (Table |3|) during this time 
again implies that either the luminosities of the two com- 
ponents are correlated, or that the underlying spectral 
shape is a more complex "single component" that hap- 
pens to mimic a power-law and a blackbody. 

The spin-down rate of IE 1048.1-5937 obtained from 
ROSAT and ASCA observations (Mereghetti 1995; Corbet 
& Mihara 1997), showed an 80% increase with respect to 
the value of — 5 x 10 -4 s yr _1 measured before 1988 with 
Ginga and EXOSAT. This increasing trend was further 
extended by the RossiXTE observation of 1996 July (P 
= 6.449769 + 0.000004 s, Mereghetti et al. 1997), yielding 
a P of 13x 10~ 4 s yr -1 after the ASCA measurement. The 
pulse period obtained with BeppoSAX lies significantly be- 
low the linear extrapolation from the previous two mea- 
surements, indicating, for the first time in IE 1048.1-5937, 
a decrease in the overall spin-down rate. 

Table || and Fig. || indicate that there is no clear corre- 
lation between the observed long term spin-down rate (P) 
and the 2-10 keV source flux. At first sight this argues 
against an accretion hypothesis, but since P is a long- 
term average and the flux an instantaneous measurement 
this comparison may not be valid. A much better compar- 
ison would be between the average flux and P during the 
same time interval. Unfortunately, due to the faintness of 
the source, no such comparison is available. In view of the 
uncertainties in the average flux a good measure of the 
time variability of the source on timescales of months to 
years would be useful. 
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